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Enhancement of bulk photovoltaic effect in topological insulators 


Liang Z. Tan and Andrew M. Rappe 
The Makineni Theoretical Laboratories, Department of Chemistry, 

University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA 

We investigate the shift current bulk photovoltaic response of materials close to a band inversion 
topological phase transition. We find that the bulk photocurrent reverses direction across the band 
inversion transition, and that its magnitude is enhanced in the vicinity of the phase transition. 
These results are demonstrated with first principles DFT calculations of BiTel and CsPbl 3 under 
hydrostatic pressure, and explained with an analytical model, suggesting that this phenomenon 
remains robust across disparate material systems. 


The field of topological insulators 0 i has proven 
to be a fruitful avenue of research in condensed matter 
physics in recent years. The symmetry protected surface 
states of topological insulators M give rise to novel and 
useful phenomena such as dissipationless transport 0 - 0 . 
Even though band topology is a ground state property, 
defined on the occupied valence bands of an insulator [8J, 
the excited state properties should be affected by topo¬ 
logical considerations as well. There is a growing body of 
research investigating the optical properties of topologi¬ 
cal insulators 0] using photoemission spectroscopy 0J3, 
Raman spectroscopy 111) 1 .’ll. nonlinear optics |14l - ll9j . 
and photocurrent generation [20fl22| . 

While the presence or absence of topological surface 
states has thus far being the main experimental route for 
detecting topological phase transitions, in this work we 
propose an optical, non-surface approach for detecting 
topological phase transitions in the bulk. This is a direct 
approach which does not rely on the quality or ease of 
detection of surface states. We consider the influence of 
band topology on the bulk photovoltaic effect (BPVE, 
also known as the photogalvanic effect) [23|, [24j], which is 
the generation of photocurrents in the bulk of a single¬ 
phase material. 


The bulk photovoltaic effect is a second-order optical 
effect which gives current densities, J r = a rs tE s E t , as a 
quadratic function of the applied electric fields. This dic¬ 
tates that the BPVE can only be observed in materials 
with broken inversion symmetry. Measurements of ap¬ 
preciable BPVE have been reported for many ferroelec¬ 
tric materials [25l428j . In the prototypical ferroelectric 
BaTiO.-. experiments HM and first-principles calcu¬ 
lations [31] show that the primary mechanism for BPVE 
is the shift current (32|. In this mechanism, carriers are 
excited into a current-carrying superposition of excited 
states. Due to its dominant effect in the bulk photo¬ 
voltaic response, we will focus on the shift current in this 
paper. 


Because optical excitations probe both the valence and 
conduction band wavefunctions, we expect the band in¬ 
version process, which is the interchange of the conduc¬ 
tion and valence band characters, to result in dramatic 
changes in the finite-frequency response functions of the 
system. We show that band inversion reverses the direc¬ 
tion of the shift current. The magnitude of the shift cur- 
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FIG. 1. Structure of BiTel displayed in top view (a), and 
side view (b). Grey: Bi, purple: Te, orange: I. Panels (c) - 
(e): Evolution of band structure across the topological phase 
transition, from the normal insulator phase (lattice constant 
a = 4.42 A) (c), to the band inversion transition (a = 4.18 A) 
(d), to the topological insulator phase (a = 4.05 A) (e). Colors 
represent the projection of the wavefunctions to Bi 6 p orbitals. 
The band character changes rapidly at the band inversion 
transition (d). 


rent is enhanced in the vicinity of the band inversion tran¬ 
sition. We illustrate these results using first-principles 
density functional theory (DFT) calculations of the po¬ 
lar topological insulators BiTel and CsPbE. Under hy¬ 
drostatic pressure, these materials undergo band inver¬ 
sions and enter topologically non-trivial phases, where 
the magnitude of shift current is enhanced to an order of 
magnitude larger than that of BaTiCU. 

Being a second-order response function, the shift cur¬ 
rent bulk photovoltaic response function includes virtual 
transitions of electrons to all unoccupied bands. 
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Unlike linear response (e.g., Fermi’s golden rnle), inter¬ 
changing the conduction band and valence band in Eq. |T| 
changes the value of the response function. In Eq. [TJ the 
sum over states includes all conduction and valence bands 
and corresponding integrals over the Brillouin zone. The 
components of the momentum operator are denoted by 
p r . The shift current can be re-expressed as a product of 
the absorption rate and the shift vector R. 


<J r st(oj) = 7re ) £( c br-|^)Hp s |c} 

cvk 

<5(w c -uj v - u)Rt(c,v,k) 

( 2 ) 

d 

Rt(c,v,k) = - — arg(c\p r \v) - [Xvt(k) - Xct{k)} (3) 

Here, x are Berry connections. The shift vector R can 
be understood as a generalized gauge invariant fc-space 
derivative of the p operator [32, [33 , and it is odd under 
the interchange of c and v bands. 

We apply this formalism to perform first-principles 
DFT calculations [Til [3?f of the shift current BPVE of 
the polar layered material BiTel. This material is a nor¬ 
mal insulator under ambient pressure, and DFT predicts 
that it becomes a topological insulator at pressures higher 
than 1.7 GPa [13 • It lacks inversion symmetry on either 
side of the topological phase transition. The strong spin- 
orbit interaction, together with the breaking of inversion 
symmetry, result in a Rashba-type splitting [39l44l| of 
the conduction and valence bands in both topological in¬ 
sulator (TI) and normal insulator (NI) phases. At the 
topological phase transition, the band gap closes between 
the lowest energy Rashba conduction band and the high¬ 
est energy Rashba valence band, in the vicinity of the A 
point in the Brillouin zone (Fig. [T]). Due to trigonal warp¬ 
ing effects, the band gap closes at six symmetry-related 
points along the A — H directions. Nevertheless, because 
the magnitude of trigonal warping is small, there is an 
approximate ring-like degeneracy of the conduction and 
valence bands at the topological phase transition at an 
energy scale above 10 meV. In general, for topological 
phase transitions in three dimensions, it is expected that 
the band gap closes over a range of pressures [13 ■ How¬ 
ever, this pressure range is very small in BiTel, appearing 
only at higher orders in k (38[, and does not affect the 


nature of our results presented below. Electronic correla¬ 
tion effects beyond DFT are likely change the transition 
pressure, but not the presence of a band inversion tran¬ 
sition [ 43 . 

We have calculated band energies and electronic 
wavefunctions within DFT, using the PBE exchange- 
correlation energy functional and norm conserving pseu¬ 
dopotentials. A 8x8x8 Monkhorst-Pack fc-grid and 
plane wave energy cutoffs of 50 Ry were used for self- 
consistent evaluation of the charge densities. A re¬ 
fined fc-point sampling near the A points equivalent to 
a 96x96x96 mesh revealed changes in the shift current 
spectrum as the system is tuned across the topological 
phase transition. For simplicity, we restrict our attention 
to the a zzz tensor component of the bulk photovoltaic re¬ 
sponse in this paper. 

We focus on the low-energy shift current spectrum, 
with contributions from the lowest-energy conduction 
and valence bands. The dominant spectral feature in 
this energy range is a peak located at the band gap en¬ 
ergy (Fig. [2j). Starting in the NI phase, an increase in 
pressure reduces the band gap, driving the shift current 
peak to lower energies while simultaneously increasing its 
magnitude. At the critical pressure for band inversion, 
the shift current peak abruptly changes sign. Further 
increasing the pressure increases the band gap, shifting 
the peak to higher energies and reducing its magnitude. 
Comparing the magnitude of the shift current peak in TI 
and NI phases, we have found that the peaks in the TI 
phase are consistently larger in magnitude than in the NI 
phase (Fig. [3]). The peaks in either phase are sharpened 
as the band inversion transition is approached. 

In brief, the change in direction of the shift current at 
the band inversion transition is a result of the shift vec¬ 
tor R flipping sign upon interchange of conduction and 
valence bands. The enhancement of the shift current is 
a consequence of the rapidly changing band characters 
at the band inversion transition (see SI-II), as measured 
by the generalized fc-derivative in the shift vector Eq. [3] 
In addition, the smaller bond lengths in the TI phase in¬ 
crease the magnitude of the momentum matrix elements, 
leading to the enhancement of the shift current in the TI 
phase over the NI phase. 

To further understand these results, we develop a low 
energy effective model describing the band inversion pro¬ 
cess. We start with a 4-band Dirac model Hq which has 
been used to describe topological phase transitions in 
centrosymmetric systems, such as semiconductor quan¬ 
tum wells and perovskites [ffij . 


Ho(k) = AI(k)r z (g> do — Rivft x ® (fc • a) (4) 

Here, t and a are pauli matrices representing orbital and 
spin degrees of freedom respectively. The second term in 
Eq. [4] describes the band structure at the point of band 
inversion, where it is gapless, and has linearly dispersing 
bands with band velocity vf ■ The first term introduces 
a band gap away from the band inversion, controlled by 
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FIG. 2. Shift current bulk photovoltaic response of BiTel (a) 
and CsPbI 3 (b). In each panel, spectra are given for a range 
of lattice constants, specified in A. 


the mass term M(k) = m — Ak 2 . We have taken the 
time-reversal operator to be T = t z ®Kia y . As inversion 
symmetry breaking is a prerequisite for shift current, we 
add an inversion symmetry breaking term Hisb , which 
preserves the time-reversal symmetry of Hq. To lowest 
order in k , all such inversion symmetry breaking terms 
can be written as t x ® cq, which introduces a polar axis 
in direction i. Without loss of generality, we choose H = 
Ho + Hisb, H isb = ht x <g> a z . H ISB is responsible for 
the Rashba split conduction and valence bands (Fig. [T]). 

Only the lowest-energy Rashba conduction and valence 
bands participate in the band inversion transition. As 
such, we reduce H to a 2 x 2 Hamiltonian H' for these 
bands by means of a canonical transformation (SI). The 
reduced Hamiltonian is 

H' = W ] HW 

M(k) A \ (5) 

A —M(k) J 


FIG. 3. Spectral parameters of the shift current response for 
BiTel. Absolute value of the peak shift current is shown in (a) 
as a function of different band gap energies, corresponding to 
different applied pressures and lattice constants. Full width 
at half maximum of the shift current peak is shown in (b) 
as a function of different band gap energies, corresponding to 
different lattice constants. Dashed lines show the prediction 
of the model developed in Eq. Q using dipole matrix elements 
fitted to the normal insulator (NI) phase. 

with band energies 

E ± {k)=± sj M{k) 2 + A 2 (6) 


where A = \J(hvFk z ) 2 + (/x — hvpk_ l) 2 . From Eq. E3 we 
see that the band inversion transition happens when the 
parameter m is tuned to m* = A[r/Yrv 2 F . The band 
gap closes in the ring k z = 0,fcj_ = /j,/Hvf , with the 
size of the ring depending on the amount of inversion 
symmetry breaking (fi). From the relative signs of the 
mass term in Eq.[I]at k « 0 and k —> oo, it can be shown 
that m > m* corresponds to the inverted band structure, 
while m < m* corresponds to the normal, uninverted 
phase 0, when A > 0. In BiTel, m can be controlled 
by applying pressure. 
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We now proceed to obtain the a zzz tensor component 
of the shift current for the effective Hamiltonian H'. In 
the basis of H , the momentum operator p is taken to be 
T x Oi , by the same argument for the form of Hjsb- The 
shift current contribution from the low energy bands is 
obtained by reducing to the basis of H’ using the trans¬ 
formation W^pW, and applying Eq. [ 2 ] (SI). 

/ p \ 2 M 

CT(w) = 7re wd \p? DOS( W )0( W -2? fl ) — 

+ (^) 

= ” e (kk ) 2 |p|2 D0S <")«<" - 

Here \p\ 2 is the dipole matrix element, assumed 
to be constant in this range of k, and M(u) = 
M(k )\ u =e (k)-E-(k)- Eq. [3 shows that the shift cur¬ 
rent lineshape is a peak located at the band gap energy 
a > = Eg. The direction of shift current is dictated by the 
sign of M . Across the band inversion transition, the mass 
term M changes sign, as does the shift current response, 
in agreement with the results of DFT (Fig. [2]). We note 
that the direction of the shift current in a real material is 
a consequence of many factors [H|; band inversion flips 
the sign of the shift current, but does not fix a particular 
direction for the current. 

The enhancement of the DFT shift current in the vicin¬ 
ity of the band inversion transition is well explained by 
this model (Fig. O, and is consistent with the correla¬ 
tion between band gap and shift current reported in [45j. 
Likewise, the sharpening of the shift current peak near 
the band inversion transition can be understood by writ¬ 
ing M/E+ = M/y/M 2 + A 2 in Eq. 0 When the mass 
term M is small, it determines the width of the peak, as 
seen from the agreement between the model and DFT in 
Fig. [3] As M is increased, it reaches a point where it be¬ 
comes comparable to the Rashba splitting energy scale, 
which limits the range of applicability of Eqs. [5]0 

We expect the qualitative nature of these results to 
remain unchanged in the presence of higher-order band 
warping effects. To demonstrate this, we calculate the 
shift current response of the halide perovskite CsPbl 3 us¬ 
ing DFT. Under hydrostatic pressure. CsPbl 3 becomes a 
ferroelectric topological insulator [4g|. The band struc¬ 
ture in the vicinity of the topological phase transition 
is similar to that of BiTel, except for a warping of the 
low energy bands, which breaks the ring-like minima 
of the Rashba bands into four shallow, degenerate val¬ 
leys (Fig. |3). Nevertheless, the shift current response of 
CSPM 3 shows similar trends as that of BiTel (Fig. 0). 
The main differences arise from the different density 
of states of BiTel (DOS(w) ~ y/u — E g ) and CsPbI 3 
(DOS(w) ~ (w — Eg) 3 ’/ 2 ), which give rise to somewhat 
different line shapes, and a weaker enhancement of shift 
current for CsPbl 3 - 

Shift current materials have been suggested as candi¬ 
dates for light-harvesting energy applications, due to ac¬ 
tive photocurrent generation in the entire bulk of the ma¬ 
terial and their above-band gap photovoltages [47^ . The 
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FIG. 4. Low energy band structure of BiTel (a) and CsPblg 
(b) near the band inversion transitions. Energy difference 
between the lowest conduction band and highest valence band 
is shown at a = 7.90 A for BiTel (a) and at a = 6.075 A for 
CsPbI 3 (b). 

BPVE is a ballistic hot carrier effect. As such, a material 
with a small band gap can still have carriers emerge with 
more than the band gap worth of energy, as they do not 
thermalize before leaving the photovoltaic material. Re¬ 
cent works [37l[48l450l| have studied how to achieve a large 
shift current magnitude in real materials, by considering 
structural modifications such as domain walls and chemi¬ 
cal alloying. Our result suggests a new way to design shift 
current materials based on their electronic topology: to 
look for materials close to a band inversion transition. 
Furthermore, the shift current enhancement reported in 
this paper may find applications in high-resolution sen¬ 
sors and switches, due to its robustness and the abrupt¬ 
ness of the change in current direction. For instance, the 
direction of the photocurrent can be controlled by small 
changes in the lattice constant, flipping directions as the 
material undergoes a band inversion transition. Complex 
functionality coupling electrical, optical, and mechanical 
degrees of freedom may be implemented simply in a bulk 
crystal of these topological shift current materials. 

We have demonstrated clear signatures of band inver¬ 
sion in the shift current response of noncentrosymmetric 
materials. As a material is tuned across a topological 
phase transition, we predict a reversal of the photocur¬ 
rent direction, sharpening of the shift current peak, and 
enhancement of the photocurrent during band inversion 
and in the inverted phase. These features are robust 
across various material systems such as the perovskite 
CsPbl 3 and the layered BiTel. Besides shedding light on 
the excited state properties of topological insulators and 
providing a bulk optical signal of band inversion, these re¬ 
sults suggest a route for enhancing the bulk photovoltaic 
effect by designing materials close to band inversion, with 
possible applications in light-harvesting devices. 
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